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Abstract. A canonical-basis formulation of the time-dependent Hartree-Fock-Bogoliubov 
(TDHFB) theory is obtained with an approximation that the pair potential is assumed to be 
diagonal in the time-dependent canonical basis. The canonical-basis formulation significantly 
reduces the computational cost. We apply the method to linear-response calculations for even- 
even nuclei. El strength distributions for proton-rich Mg isotopes are systematically calculated. 
The calculation suggests strong Landau damping of giant dipole resonance for drip-line nuclei. 



1. Introduction 

The time-dependent Hartree-Fock (TDHF) theory has been extensively utilized to study nuclear 
many-body dynamics [Ij. Recently, it has been revisited with modern energy density functionals 
and more accurate description of nuclear properties has been achieved [2J O HI [SJ El E] . The 
TDHF theory uses only occupied orbitals, number of which is equal to the number of particles 
(N). However, it neglects the residual interactions in particle-particle and hole-hole channels, 
which are important for properties of open-shell heavy nuclei. It is well-known that the time- 
dependent Hartree-Fock-Bogoliubov (TDHFB) theory [8j properly takes into account the pairing 
correlations. The TDHFB equation is formulated in a similar manner to the TDHF, however 
it requires us to calculate the time evolution of quasi-particle orbitals, number of which is, in 
principle, infinite. Therefore, the practical calculations with the TDHFB are very limited [9lll0|. 

Very recently, we have proposed a possible approximation for the TDHFB theory This is 
a time-dependent version of the BCS approximation [12J for the Hartree-Fock-Bogoliubov theory. 
Namely, we neglect off-diagonal elements of the pair potential in the time-dependent canonical 
basis. We show that this approximation results in significant reduction of the computational 
task. We call the equations obtained with this approximation, "Canonical-basis TDHFB" (Cb- 
TDHFB) equations. We apply the method to the linear-response calculations using the full 
Skyrme functionals of the parameter set of SkM*, and discuss properties of the El strength 
distribution in even-even Mg isotopes. 

The paper is organized as follows. In Sec. [2j we present the basic equations of the present 
method and their derivation. It is emphasized that the basic equations possess a gauge 
invariance. A schematic choice of the pairing functional leads to violation of the gauge invariance, 
which requires us to choose a specific gauge to minimize the violation. In Sec. El we present 
numerical results of the real-time calculations of the linear response for Mg isotopes. Finally, 
the conclusion is given in Sec. 01 



2. Formalism of the Cb-TDHFB theory 

In this section, we show the basic equations of the Cb-TDHFB and recapitulate their derivation. 
See Ref. [11] for more details. 



2.1. Basic equations 

Let us first show the Cb-TDHFB equations we derive in the followings. 

iftlMt)) = (M*) - vk(t))\Mt)), ig- t \fa(t)) = (h(t) - v- k (t))\h(t)), (la) 

ij t Pk(t) = Kk(t)A%(t) - K* k (t)A k (t), (lb) 
ij t K k(t) = (Vk(t) + Vk(t)) K k (t) + A fe (t) (2p k (t) - 1) . (lc) 

These basic equations determine the time evolution of the canonical states, \4>k(t)) and |<^(i)}, 
their occupation, pk(t), and pair probabilities, Kk(t). The real functions of time, r]k(t) and nj,{t), 
are arbitrary and associated with the gauge degrees of freedom. The time-dependent pairing 
gaps, Afc(i), which are given in Eq. (|12p . are similar to the BCS pairing gap |12j except for 
the fact that the canonical pair of states are no longer related to each other by time reversal. 
Although we use the same symbols, (/?, k, A), for matrixes in Eqs. (J3aj) and ([3b]), the quantities in 
the Cb-TDHFB equations are only their diagonal elements with a single index for the canonical 
states k. It should be noted that similar equations can be found in Ref. [T3] for a simple pairing 
energy functional. 

2.2. Derivation of the basic equations 

We now derive the Cb-TDHFB equations starting from the generalized density-matrix formalism. 
The TDHFB equation can be written in terms of the generalized density matrix R(t) as [8j 

i^R(t) = [H(t),R(t)]. (2) 

This is equivalent to the following equations for one-body density matrix, p(t), and the pairing- 
tensor matrix, n(t). 

i^p(t) = [h(t), P (t)] + K(t)A*(t)-A(t)K*(t), (3a) 
*^«(t) = h(t)K{t) + K{t)h*(t) + A{t){l-p*{t))-p{t)A{t). (3b) 

Here, hit) and A(t) are single-particle Hamiltonian and pair potential, respectively. 

At each instant of time, we may diagonalize the density operator p in the orthonormal 
canonical basis, {4>k(t), <t>k(t)} with the occupation probabilities pk- Then, the TDHFB state is 
expressed in the canonical (BCS) form as 

i*(*)> = n + ^(*)4w#)} i°>- ( 4 ) 

fc>0 

For the canonical states, we use the alphabetic indexes such as k for half of the total space 
indicated by k > 0. For each state with k > 0, there exists a "paired" state k < which 
is orthogonal to all the states with k > 0. The set of states {(/>&, <?^} generate the whole 



single-particle space. We use the Greek letters p, u, ■ ■ ■ for indexes of an adopted representation 
(complete set) for the single-particle states. Using the following notations, 

«H0 fc (t)0ifc(t)» = (v\Mt))(v\h(t))-(v\Mt))(v\Mt)), (5) 
= \Mt))(Mt)\ + \^k(t))(Mt)\, (6) 

the density and the pairing-tensor matrixes are expressed as 

P»u(t) = J2p k (t){fi\n k (t)\u), (7) 

fc>0 

K^(t) = (8) 



fc>0 

|2 



where p k (t) = \v k (t)\ and K k (t) = u* k (t)v k (t). It should be noted that the canonical pair of 
states, \4> k {t)) and \<t> k (t)}, are assumed to be orthonormal but not necessarily related with each 
other by the time reversal, \<p k ) / T\<p k ). 

We can invert Eqs. (|7|) and (JSj) for p k and K k , 

Pkit) = ^(^(t)|/i)^(t)(H^(t)} = ^(^WWv(t)(^^)), (9) 

fit/ flU 

«*(*) = lj2^ k ^(t)\HKu(t). (io) 
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The derivative of p k (t) with respect to time t leads to 

flU 



^{ Kfe (i)A^(t){(^|0 fe (t)^(t)))+4(t)A^(t){(^(i)^(i)|H)}. (11) 
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We used the assumption of norm conservation and the TDHFB equation ()3a|) . This can be 
rewritten in the simple form of Eq. (jlb[) with the definition of the pairing gap, 

A k(t) = -^X) A ^(*)«^(*)^(*)iM>- ( 12 ) 

In the same way, we evaluate the time derivative of K k (t) as 

4 Kfe( * } ^E^^^^wi^^+^w O^r |0fe(t)> + 0^ (t)> ) • (13) 

Then, using the TDHFB equation (|3bp . we obtain Eq. (fTc|) with the real gauge functions 



m (t) EE <^)|fc(t)|0 fc (t)) + i(-^\ct, k (t)). (14) 

These functions control time dependence of phase for the canonical states, which are basically 
arbitrary. 



So far, the derivation is based on the TDHFB equations, and no approximation beyond the 
TDHFB is introduced. However, to obtain simple equations for time evolution of the canonical 
basis, we need to introduce an assumption (approximation) that the pair potential is written as 

fe>0 

This satisfies Eq. (f!2j) . but in general, Eq. ([12]) can not be inverted because the two-particle 
states {(pkfik) do not span the whole space. In other words, we only take into account the 
pair potential of the "diagonal" parts in the canonical basis, A k j = —A k 5 k i. In the stationary 
limit (\4>k) = T\4>k)), this is equivalent to the ordinary BCS approximation |12| . With the 
approximation of Eq. (fT5j) . it is easy to see that the TDHFB equations, (f3a|) and (|3bj) . are 
consistent with Eqs. ([TJ. 

2. 3. Properties of the Cb- TDHFB equations 

The Cb-TDHFB equations, ([1]), are invariant with respect to the gauge transformation with 
arbitrary real functions, 9 k {t) and k (t). 

\4> k ) e»*®\<f> k ) and \cj> k ) -+ e**®^) (16) 
Kk ^ e -i(0k(t)+e- k (t)) Kk and A k -y e -^(*)+%W) A k (17) 

simultaneously with 

Vk(t)^Vk(t) + ^- and Vk(t)^ Tk (t) + ^. 

The phase relations of Eq. (|17p are obtained from Eqs. (jlpp and ()12p . It is now clear that the 
arbitrary real functions, n k (t) and rj k (t), control time evolution of the phases of \4> k (t)), \4> k (t)), 
Kk(t), and A k (t). 

In addition to the gauge invariance, the Cb-TDHFB equations possess the following 
properties. 

(i) Conservation law 

(a) Conservation of orthonormal property of the canonical states 

(b) Conservation of average particle number 

(c) Conservation of average total energy 

(ii) The stationary solution corresponds to the HF+BCS solution. 

(iii) Small-amplitude limit 

(a) The Nambu-Goldstone modes are zero-energy normal-mode solutions. 

(b) If the ground state is in the normal phase, the equations are identical to the particle- 
hole, particle-particle, and hole-hole RPA with the BCS approximation. 

2.4- Energy Junctionals and numerical procedure 

We adopt a Skyrme functional with the SkM* parameter set for the particle- hole channels. For 
the pairing energy functional, we adopt a simple functional of a form 

E S) = " E G klK * k (t) Kl (t),= -Y,<(t)A k (t), A k (t) =Y,G k m(t), (18) 

k,l>0 k>0 l>0 

where G k i = G f (e k ) f (e®) with G = 0.6 MeV. The smooth cut-off function f(e k ), whose explicit 
form can be found in Ref. [11J . depends on the single-particle energy of the canonical state k at 



Table 1. Ground-state properties of Mg isotopes calculated with the SkM* functional; 
quadrupole deformation parameters (/3,7), pairing gaps for neutrons and protons (A n ,A p ), 
chemical potentials for neutrons and protons (A n , A p ). In the case of normal phase (A = 0), we 
define the chemical potential as the single-particle energy of the highest occupied orbital. The 
pairing gaps and chemical potentials are given in units of MeV. 





P 


7 


A n 


A p 


— A n 


— A p 


18 Mg 


0.31 


0° 


0.0 


0.0 


25.59 


0.20 


20 Mg 


0.0 




0.0 


1.13 


20.53 


2.83 


22 Mg 


0.38 


0° 


0.0 


0.0 


16.31 


6.42 


24 Mg 


0.39 


0° 


0.0 


0.0 


14.12 


9.51 


26 Mg 


0.20 


54° 


0.0 


0.86 


13.08 


11.23 


28 M g 


0.0 




0.0 


1.03 


9.21 


13.30 



the HF+BCS ground state. A drawback of the simple functional (|18p is that it violates the gauge 
invariance, which results in breakdown of some of nice properties shown in Sec. 12.31 However, 
it is shown that all these properties can be recovered by choosing a special gauge condition |11] 

vk(t) = e k (t) = (Mt)\h(t)\Mt)), m(t) = e- k (t) = (h(t)\Ht)\h(t)). (19) 

For numerical calculations, we extended the computer program of the TDHF in the three- 
dimensional coordinate-space representation [2] to include the pairing correlations. The ground 
state is first constructed by the HF+BCS calculation. Then, we solve the Cb-TDHFB equations 
in real time, under a weak impulse isovector dipole field, yielding a time-dependent El moment, 
Dgiit). To obtain the El strength distribution, we perform the spectral analysis with an 
exponential smoothing with V = 1 MeV: DEi(t) — > DEi(t)e~ rt ^ 2 . The readers should refer to 
Ref. |11| for more details. 

3. Electric dipole strength distribution in proton-rich Mg isotopes 

In Table [lj the ground-state deformations, pairing gaps, and chemical potentials are listed for 
stable to proton-rich Mg isotopes. In the present calculation with SkM*, 18 Mg turns out to be 
bound with a small binding energy of 200 keV. Thus, we include this nucleus in our calculation 
as a "fictitious" proton halo nucleus. The neutron pairing gap is absent for all these nuclei. The 
proton gap also vanishes for nuclei with prolate shapes, (3 = 0.3 ~ 0.4. 

Calculated El strength distributions are shown in Fig. [TJ The double-peak structure of the 
giant dipole resonance (GDR) due to the deformation splitting can be seen in 22 ' 24 Mg. The 
K = peak is located around 15 MeV and the K = 1 is near 22 MeV. In contrast, for 18 Mg, 
although the ground state is deformed in a prolate shape with (3 « 0.3, the double-peak structure 
is not clearly seen. In this nucleus, the El strength of both the K = and K = 1 components are 
fragmented into a wide range of energy. Previously, we calculated the El strength distribution 
in neutron-rich Mg isotopes |llj . and found the similar damping effects near the drip line. This 
strong Landau damping near the drip lines may be understood by the high level density of 
one-particle-one-hole (lplh) states near the GDR energy. 

The low-energy El strength in E < 10 MeV is negligible for stable nuclei ( 24,26 Mg). For the 
neutron-rich side, we see a small low-energy peak in 28 Mg. Since the neutron separation energy 
is still sizable (about 9 MeV), we assume that this is due to the occupation of the neutron si/ 2 
orbital which is spatially extended. For the proton-rich (neutron-deficient) side, the low-energy 
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Figure 1. Calculated El strength distribution for even-even Mg isotopes (JV = 8 ~ 16). For 
deformed nuclei, the total strength is decomposed into K = (green dashed line) and K = 1 
(blue dotted line) components. The z-axis is chosen as the symmetry axis for axially deformed 
cases. The smoothing parameter of T = 1 MeV is used. 



strength is seen in 18,20 Mg. This should be due to the weak binding of the last-occupied proton 
d 5 / 2 orbitals, since the calculated proton separation energies are less than 3 MeV for these nuclei. 

4. Summary 

We presented an approximate and feasible approach to the TDHFB; canonical-basis TDHFB 
method. Since the number of the canonical states we need to calculate is the same order as the 
particle number, this method significantly reduces the computational task of the TDHFB. We 
calculated the El strength distribution in proton-rich Mg isotopes, using the real-time real-space 
method. We found a strong Landau damping effect in the drip-line nuclei, that may be related 
to high level density of the background lplh states with negative parity. The calculation also 
indicates an increase of the low-energy El strength as the nucleus approaches the proton drip 
line. 
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